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ABSTRACT 

Frequency parameter plane techniques are developed for 
the analysis and synthesis of both linear and nonlinear 
systems. The transmission function of a dynamic system 
contains two parameters which are adjusted to provide 
desired performance. These parameters represent elements 
of the system. 

The four dimensional paremeter space relating magnitude, 
frequency, and the two parameters may be represented by seven 
two dimensional projections of the parametric curves. 

Computer programs are developed for the computation, and 
graphs of results are obtained. Interpretations of the graphs 
can be used for design of systems such as tank stabilizers. 
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I. INTRODUCTION 



The analysis and synthesis of control systems may be 
accomplished using the frequency parameter plane technique. 

A primary consideration in control system design is the 
adjustment of the system frequency response. To achieve 
this, the frequency parameter plane technique may be employed 
and is based on a graphical procedure for considering four 
parameters, alpha, beta, frequency and magnitude of the 
transmission functions. The objective is to choose values 
for alpha and beta (or functional relationships if the design 
is to be nonlinear) to provide the desired frequency response 
characteristics. This choice is possible through proper use 
of the frequency parameter plane curve families. 

The goal of this thesis is design using analysis techniques. 
The parameter frequency response method is used to synthesis 
and design systems. A steam-engine governor is used in the 
analysis as an illustration. The antiroll tank stabilizer of 
a ship is also designed using the frequency parameter plane, 
for both passive and active tanks. The tank size and active 
device characteristics are represented by the parameters, 
which are varied. Parameter plane curve families are used 
in the design procedure presented in later chapters. 
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II . PARAMETER PLANE ANALYSIS 



Dynamic systems may be either linear or nonlinear. 

Parameter plane analysis permits us to study dynamic behavior 
by treating system parameters as variables. One first 
determines the characteristic polynomial of the system which 
may be in the form 

N b M M_1 

F(s) = E aj^ s'" = a^jS " ' + + a^S+a^ = 0 

k * 0 

( 2 . 1 ) 

where a|^= coefficients which are functions of the parameters. 

S = o-+joj’ = _^2 = w^e'^® 

By substituting the value of S into equation 2.1, 
equation 2.2 is obtained. 

F(s) = E a. 0)1" e^*^® (2.2) 

k=0 " 

Since equation 2.2 has complex coefficients it can be converted 
into two equations by requiring reals and imaginaries to be 
zero independently. 

Note that S*^ = (Cos k6 + j Sin k0) 

= + j 

where = oj|^Cos k0 = oj^Cos(k Cos"^?) = { -1 ) ( ?) 

Yk = w|^Sin k0 = w|^Sin(k Cos'^c) = ( -1 ) '^aj^U|^( ?) 
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T|^(^') and are Chebyshev polynomials with values 

as shown in Table 1, 2, and 3. 

It follows that 

N k k 

2 a. [(-D^con. (c)] = 0 
k=0 

^ k k - 

2 a. [(-D'^w'^U. (c)] = 0 (2.3) 

k=0 ~ 



are simultaneous equations which may be solved for two 
undefined parameters in the aj^. 

For many practical problems the algebraic form of the 
coefficients is one of the following two cases: 

1 . a|^ = ^k*^ ^k 

2. aj^ = bj^a + Cj^3 + hj^a3 + dj^ 

Case 1 can be solved by Cramer's rule and the results 
for a and 3 are 



a = 



3 = 



^ 1°2 ~ ^ 2°1 
^ 1^2 ■ ^ 2^1 

^ 2 °! - ^ 2 h 
^ 1^2 ■ ^ 2^1 



where 



N 



N 



B, V=0<-’) 



N 



N 



Cl Mo'-') Vn Tk(C); Cj 
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Case 2 by substituting equation 2.3 rearranges into 
the form. 

aB-j + 3C-| + a3H-| + D-j = 0 

082 + 3 C 2 + a3H2 + D 2 = 0 (2.4) 

where the definitions of case 1 still apply, and 

Eliminating 3 the solutions are obtained 

_ -e±/e^- 4ac , ^ 

“ 1,2 2a (2.5a) 

and corresponding values for 3 are 

B,a, o+D, E[„a, o+D, 

a = ‘ ^ ^ ' (o 

1,2 ^i“i,2'^^l ^2“l,2'^^2 

Eliminating a, the following solution for 3 



and the corresponding values for a are 

“1.2 - - H, 6 , ^ 3.81 - - 836 ,^ 2^83 (2.6b) 



where the notation is the following 



As can be seen, there are generally two pairs of solu- 
tions (a, 3) in terms of and C. The solutions are given 

by equations 2.5a and 2.5b provided a 0 . If a = 0 the 
solutions are given by equations 2.6a and 2.6b provided 
b 0. When a = b = 0 case 2 reduces to case 1. 

11 1 ustrati ons I . A feedback control system is as shown 





C 




The closed-loop transfer function is: 



K 



T(s) = 



+ 3S^ + (2+K^K2)S + K 



The characteristic equation of the system is 



+ 3S^ + (2+K^K2)S + = 0 



Let a = K^K 2 

6 = K, 



The equation becomes: 



+ 3S^ + (2+d)S +6=0 



From the parameter plane method 






C 



2 



0 
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= (4?^-3?)(-o:^^) + (2 c2.i)(3o)^2j ^ (.2?u)J 

^2 (4?^-l ) ( -to^^) + (2c)(3Wj^^) + (-2 o3^) 

SO, the solutions for a and 3 are 
a =-4w^^C^ + 6co^c + - 2 

B = 

The c-curves and and u^-curves on the a-3 plane are shown 
on Fi gure 2.1. 

I 1 1 ustrati on II . The system is the same as illustration I 
let = 3 

K 2 a 

The characteristic equation becomes: 

+ 3S^ + (2+a3)S +3=0 

Then 

Bj = 0 B 2 = 0 

= 1 C2 = 0 

Hi = -Coj^ H2 = 

= (3?-4.C^,) + (6C--3)o3^^ " ^ 

D 2 = (1-4 + 6Co 3^^ - 2o)^ 

And the other definitions are 

a=0; b=u^; c=0; 

d = + 6 + (3w^-2w^)C - 3w^ 

e = f = 2u^ ? - 3to^ 
n n 
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So only one solution for a and g in equations 2.6a and 2.6b 
is obtained: 



The curves are shown on Figure 2.2 

Consider the curves on Figures 2.1 and 2.2: the roots of 
the characteristic equation are defined as functions of a and 
g. If we choose g = 40 and a = 20 on Figure 2.1, then 



therefore, the roots of the system are: 

S = - 0.5 ± J4.47 

and the corresponding value for the real root is r2. In like 
manner on Figure 2.2, if we desire roots at c = 0.111 and 
03^ = 4.5 we found a = 0.5 and g = 40. 

On the other hand, if ^ = 0.2 and = 5 the roots are 
-l±j4.9 and -1. From the Figure we find: 
g = = 25 

a = 25 (Figure 2.1) or 1.0 (Figure 2.2) therefore K 2 = 1 .,0 
Parameter plane methods are very convenient. The a-g 
plane can be used to find the system parameters, and/or roots 
of the characteristic equation at any defined values of the 
parameter. The stability of the system can be investigated 
at the line c=0 which corresponds to the imaginary axis on 
the S-plane. 






3 2 

2a)'"c; - 3o) 

n n 



“n ' '*•5 

? = 0.111 



16 



i 



aulc !• uMt:u^:>ncv ruriuuiuri:> i V S V » 



O 



o 

II 



o 



+ 

I — CO 
+ XJ' 

C\J o 

XJ' 1£) 

o I— 

XJ' 1 

00 I lO 
I ^ XJ' 













p— 


CO 


XJ' 


CSJ 












1 


XJ' 


o 














CSJ 


o 




LD 












XJ' 


00 


CSJ 


+ 










XJ' 




+ 


+ 










p— 


ID 


CSJ 


in 


LD 


XJ' 








+ 


+ 


+ 


kJ' 


kJ' 










CSJ 


CO 




CSJ 


00 


CSJ 








XJ' 


XJ' 


XJ' 


a> 




o 








CSJ 


CSJ 


o 


I — 




I — 






XJ' 


P-* 


CO 


00 


1 


1 


1 




p— 




1 


1 


1 




00 


cn 




1 


1 




LO 


LD 


XJ' 


XJ' 


XJ' 




CJ 


CO 


XJ' 


XJ' 


XJ' 


00 


LD 


CSJ 


XJ' 


XJ' 


XJ' 


K£> 


CSJ 




CSJ 


in 


p— 


CJ 




00 


•“ 


CO 


so 


n- 


CSJ 


in 


II 


II 


II 


II 


II 


II 


II 


II 


II 


















o 


C\J 


CO 




in 


so 




00 


cr» 


1 — 


ZD 


ZD 


ZD 


ZD 


ZD 


ZD 


ZD 


ZD 


ZD 



o 



I 

CJ 

XJ' 

o 

lO 

+ 

cr> XJ' 

















+ 


o 














+ 


CO 


o 














CSJ 


XJ' 
















kj' 


o 


1 














CSJ 


CSJ 


LD 














CO 


r— 


XJ' 














1 


1 


o 












r>- 




in 


CSJ 










r— 


1 


XJ' 


kj' 


I — 










1 


CO 


o 


CSJ 


p— 










CSJ 


kJ' 


LD 


CO 


+ 










XJ' 


ID 


p— 




00 








XJ' 


00 


in 


+ 


+ 


XJ' 








m 


p— 


■1- 


ID 




o 








+ 


+ 


in 


XJ' 


XJ' 


00 






p— 


CO 




XJ' 


ID 


ID 


CSJ 






+ 


XJ' 


XJ' 


CSJ 


in 




I — 






CSJ 


o 


00 


1 — 


CSJ 


in 


1 




XJ' 


XJ' 


CSJ 




I — 


1 


1 


o 


p— 


CO 


00 


1 


1 


1 


00 


cr» 


I — 


1 


1 


1 


in 


so 




XJ' 


kj' 


kJ' 


CSJ 


CO 




XJ' 


XJ' 


XJ' 


00 


LO 


CSJ 


XJ' 


XJ' 


kj' 


LD 


CSJ 




CSJ 


in 


p— 


XJ' CSJ 




00 


*— 


CO 


so 


>— 


CSJ 


m 


II II 


II 


II 


II 


II 


II 


11 


II 


II 


















o 


I— CSJ 


CO 




in 


so 




00 


<Tk 


1— 


h- 1 — 


1— 


1“ 


1— 


H- 


K- 


1— 


1— 


1— 



17 



Table 2 . Functions T |^(?) 



5 To 


Tl 


T2 


T3 


T4 


T5 


Te 


0.00 


0.00 


- 1.000 


0.0000 


1 .00000 


0.000000 


- 1.0000000 


0.05 


0.05 


- 0.995 


- 0.1495 


0.98005 


0.347505 


- 0.9452995 


0.10 


0.10 


- 0.980 


- 0.2960 


0.92080 


0.480160 


- 0.8247680 


0.15 


0.15 


- 0.955 


- 0.4365 


0.82405 


0.683715 


- 0.6189355 


0.20 


0.20 


- 0.920 


- 0.5680 


0.69280 


0.845120 


- 0.3547520 


0.25 


0.25 


- 0.875 


- 0.6875 


0.53125 


0.953125 


- 0.0546875 


0.30 


0.30 


- 0.820 


- 0.7920 


0.34480 


0.998880 


0.2545280 


0.35 


0.35 


- 0.755 


- 0.8785 


0.14005 


0.976535 


0.5435245 


0.40 


0.40 


- 0.680 


- 0.9440 


- 0.07520 


0.883840 


0.7822720 


0.45 


0.45 


- 0.595 


- 0.9855 


- 0.29195 


0.722745 


0.9424205 


0.50 1 


0.50 


- 0.500 


- 1 .0000 


- 0.50000 


0.500000 


1 .0000000 


0.55 


0.55 


- 0.395 


- 0.9845 


- 0.68795 


0.227755 


0.9384805 


0.60 


0.60 


- 0.280 


- 0.9360 


- 0.84320 


- 0.075840 


0.7521920 


0.65 


0.65 


- 0.155 


- 0.8515 


- 0.95195 


- 0.386035 


0.4501045 


0.70 


0.70 


- 0.020 


- 0.7280 


- 0.99920 


- 0.670880 


0.0599680 


0.75 


0.75 


0.125 


- 0.5625 


- 0.96875 


- 0.890625 


- 0.3671875 


0.80 


0.80 


0.280 


- 0.3520 


- 0.84320 


- 0.997120 


- 0.7521920 


0.85 


0.85 


0.445 


- 0.0935 


- 0.60395 


- 0.933215 


- 0.9825155 


0.90 


0.90 


0.620 


0.2160 


- 0.23120 


- 0.632160 


- 0.9066880 


0.95 


0.95 


0.805 


0.5795 


0.29605 


- 0.017005 


- 0.3283595 


1.00 


1.00 


1 .000 


1 .0000 


1.00000 


1.000000 


1 .0000000 


Tb 




T4 


T5 




TlO 





0.00000000 1 
- 0.44203495 0 

- 0.64511360 0 

- 0.86939565 0 

- 0.98702080 -0 
- 0.98046875 -0 
- 0.84616320 -0 
- 0.59606785 -0 
- 0.25802240 -0 
0.12543345 
0.50000000 
0.80457355 
0.97847040 
0-97117085 
0-75483520 
0-33984375 
- 0.20638720 
- 0.73706135 
- 0.99987840 
- 0.60687805 
1.00000000 



000000000 
901 096005 
695745280 
3581 1 6805 
040056320 
435546875 
762225920 
960771995 
988689920 
- 0.829530395 
- 0.500000000 
- 0.053449595 
0.421972480 
0.812407605 
0.996801280 
0.876953125 
0.421972480 
- 0.270488795 
- 0.893093120 
- 0.824708795 
1 .000000000 



0.0000000000 
0.5321 445505 
0.7842626560 
0.9768306915 
0.9709982720 
0.76269531 25 
0.3888276480 
- 0.0764725465 
- 0.5329295360 
- 0.87201 08055 
-T .0000000000 
- 0.8633668104 
- 0.4721034240 
0.0849590365 
0.6406865920 
0.9755859375 
0.8815431680 
0.2772303985 
- 0.6076892160 
- 0.9600686605 
1 .0000000000 



- 1 .00000000000 

- 1 .84788154995 
- 0- 53889274880 
- 0- 06506759755 
0-42845562880 
0-81689453125 
0- 99552250880 
0-90724121245 
0-46234629120 
0-04472067005 
■ 0- 50000000000 
- 0-98625531995 
- 0-98849658880 
- 0-74178047245 
- 0.09984005120 
0.58642578125 
0.98849658880 
0.74178047245 
- 0.20074746880 
- 0.99942165995 
1 - 00000000000 
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Table 3 . 



Functions U|^(c) 



^ u -1 


"l 


^2 


^3 


^4 


^5 


^6 


0.00 




0.0 


- 1.00 


0.000 


1.0000 


0.00000 


0.05 




0.1 


- 0.99 


- 0.199 


0.9701 


0.29601 


0.10 




0.2 


- 0.96 


- 0.392 


0.8816 


0.56832 


0.15 




0.3 


- 0.91 


- 0.573 


0.7381 


0.79443 


0.20 




0.4 


- 0.84 


- 0.736 


0.5456 


0.95424 


0.25 




0.5 


- 0.75 


- 0.875 


0.3125 


1 .03125 


0.30 




0.6 


- 0.64 


- 0.984 


- 0.0496 


1.01 376 


0.35 




0.7 


- 0.51 


- 0.057 


- 0.2299 


0.89607 


0.40 




0.8 


- 0.36 


- 0.088 


- 0.5104 


0.67968 


0.45 




0.9 


- 0.19 


- 0.071 


- 1 . 7739 


0.37449 


0.50 -1 


0 1 


1.0 


0.00 


- 1 .000 


- 1.0000 


0.00000 


0.55 




1.1 


0.21 


- 0.869 


- 1 .1659 


- 0.41 349 


0.60 
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Figure 2.1 Root Parameter Plane Curves for illustration 1. 
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Figure 2.2 Root Parameter Plane Curves for Illustration 2. 
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Ill . PARAMETER FREQUENCY RESPONSE TECHNIQUES 



A. INTRODUCTION 

Both linear and non linear dynamic systems may be sub- 
jected to periodic inputs. It is important therefore that 
we be able to analyze such systems in the frequency domain. 

If the system is linear there also exists a correlation 
between frequency response and time response as guaranteed 
by the Fourier integral which permits us to estimate the 
transient response from the frequency response. Thus the 
frequency domain may be used to design systems for non- 
periodic inputs. In this chapter frequency response analysis 
of linear systems on a frequency parameter plane is discussed 
i n detail. 



B. ANALYSIS 



A linear system can be described mathematically by a 
transmission function: 



T(s) 



N(s) 

DTs7 



where the system may be either open-loop or closed-loop and 
N(s) and D(s) are polynomials. In polynomial notation 

• k 



T(s) = 



m 

E 9 1, S 

k=0 

n 

E b„S^ 
q=0 ^ 



where, a's and b's are coefficients which are functions of 
the system parameters. 
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If we choose two systems parameters a and 3 and if we 



consider only sinusoidal signals the frequency response 

transmission function becomes 

T/ - N( ja).g,3) 

■ D( jo),a,3) 

aQ + a-|ja) + a2(jo))^+ + ajj^Cjcj)'" 

b(j + b^jo) + b2(jo))^+ + b^(jo))" 

Lap - a2w + a^o) - J+jLa-jCj-a^w ■’■agw - 

[bo - b2W^ + b^oj^ - ]+j [b-jOj-b^w^+bgOj^- 









In abbreviated notation 

Ng(w,a,3) + jNQ(o),a,3) 

" Dg(w,a,3) + jD^(aj,a,3) 

Si mi 1 arly 

Ng(w,a,3) - j NQ(w,a,3) 

T(-j.o)) - _ j DQ((jj,a,3) 

Let the magnitude be represented by M: 



M = |T(jo)) I 

= |T(ja))|2 = lT(jo))| |T(-jo))l 

therefore 

2 [Ng(w ,a,3) - [NQ(oj,a,3) 

M = P ^ 

[Dg(w,a,3)] - [DQ(o),a,3)] 



from whi ch 



M 



2 



r 2 4 

[aQ-a2W +a^w 

[bo-b20j^+b^oj^ 



+ [a^o)-a2W^+ 

+ [b^w-b2W^+ 



2m 

E 

_ k=0 
2n 
E 

q=0 



^ 9m 1 k-1 

y I K 2m-i — p— 

(-l)^a. 0 )“^ + E (-1) ^ 

(even) k = l (odd) 






(-1)2. 

(even) q “ 



2n-l q-1 

q + 2 (-1) 2 ■ q 

q=l (odd) °q “ 
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Then 



m2 = 



m 

2 (-l)''a^ 0) 

k=0 



k " 2k 
k 



l" (-1)V.2q 

q=0 ^ 



where, the ' s and ' s are functions of a and 3. So, there 
are four variables M, to, ct, and 3 in this equation. The 
equation can be plotted by fixing one variable and using one 
variable as a parameter. Seven different combinations are 
possible . 

Exampl e . Consider the system of a Steam engine Governor 
[9] for which the transmission function can be written as: 



i = 
<!>. ■ 



Vi 

I m 



ml Imlm Imimlm Im 



where t|) = instantaneous speed response. 

(J)^ = command speed. 

m = equivalent mass of the governor sleeve. 

Cg = damping coefficient at the governor sleeve. 
Kg = stiffness of governor spring. 

= increase in the upward force on the governor 
si ee ve . 

= increase in steam torque. 

I = equivalent moment of inertia of the engine. 

Cg = damping coefficient of the engine. 

Kg = stiffness of engine spring. 



24 



From studying the mechanical characteristics [10], 
possible values of the parameters of the system are sub- 
stituted in the equation. Let and varied by adjusting 

the governor spring and the valves for steam flow respectively. 
Then the equation becomes: 

. 626aS 

± = 

***0 S^+1 0.5S^+( 30 + 253 )S^+(1 37. 53+6 26a+l 2.5) S+62.5B 

where a = C 2 




then 

|y|2 _ [626aa)]^ 

[w^- ( 30+253) 01^+62. 53]^+[-10.5o)^-{l 37. 53+626a+l 2. 5)(u]^ 

From the curves of Figure 3.1 several useful bits of 
information can be obtained by inspection: 

(a) At any intersection between two constant o) curves 
the coordinates define the a-3 pair, and the values of oi of 
the intersecting curves define the bandwidth (since the mag- 
ni tude is -3 db) . 

(b) A locus of constant oi define a set of a, 3 pairs 
which are potentially useful in self adaptive applications, 
i.e., if, in a given system, 3 is subject to change, values 
of a can be determined such that a magnitude M = can be 
maintained at w = Wi , by determining the values assumed by 3 
and providing a suitable adjustment in a. 

(c) If a number of families of curves are generated with 
coordinates as in Figure 3.1 but with different M values, the 
entire frequency response can be determined for any chosen 

a, 3 pair. 



25 



From the curves of Figure 3.2a, b and c, it can be seen 

that: 

(a) Each frequency has only one a, S> pair which gives 
the maximum magnitude. 

i.e., at oj = 10, a = 1.0, B = 3.0 (Figure 3.2a) 

(t) = 12, a = 1.4, B = 4.5 (Figure 3.2b) 

It) = 14, a = 1.75, B = 6.6 (Figure 3.2c) 

(b) The locus of the maximum magnitude point is a 
straight line when frequency is increased. 

(c) The axis which is drawn through the ellipses gives 
the values of B (and/or a) and magnitude which are coordinates 
of points on the envelope of the a versus magnitude curve as 
in Figure 3.3. 

On Figure 3.3 are curves of a versus magnitude. These 
curves as drawn show the envelope of the family of curves 
that corresponds to the axis of the ellipses on Figure 3.2. 

This envelope bounds all possible magnitudes at a frequency 
of 12 rad/sec, for all B* 

Figure 3.4, which presents a family of constant magnitude 
curves on a versus frequency coordinates (B = constant) pro- 
vides a global view of the possible frequency responses of 
the system: 

(a) We can plot the Bode diagram by reading the frequency 
and magnitude at the intersections along any line of constant a. 

(b) The maximum magnitude for the given values of B is 
found at the center of the family of ellipses, and also 
increases as a straight line when B is increased. 



26 



(c) At constant B (3=5.0), the curve of magnitude 
is -3 db., it shows the bandwidth is between the intersection 
of the curve and any line of constant a, such as at a = 2 
the bandwidth is 10.5 rad/sec. 

A plot of 3 versus frequency with constant a may be used 
in the same way. 

Curves of a versus magnitude are shown on Figure 3.5 
with 3 constant and frequency the family parameter. 

(a) The Bode diagram can be plotted at any a by evaluating 
magnitude and frequency along a horizontal line. 

(b) At the intersections of two frequency curves the 

Bode diagram for that value of a will have the same magnitude 

at each of the frequencies. 

(c) Each frequency curve defines one value of a (3 is 
constant) for which the magnitude is maximum at that frequency. 

(d) If two frequency curves intersect each other at 

the magnitude -3 db., the baridwidth at that value of a is 

defined by the two frequencies. 

Figure 3.6 plots a versus frequency with 3 as the para- 
meter. Each 3-curve is intersected twice by each a line, 
thus defining two values of frequency for the one value of 
the a, 3 pair, these two values of frequency give the band- 
width of the system (since the magnitude is -3 db). 

Figure 3.7 shows a family of Bode diagrams with a as the 
family parameter and 3 constant. It may be used to select a 
value of a, or to define a desired variation in a as a function 
of frequency. 
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Each curve family is seen to contain useful information. 
Depnding on the type of information needed some curve families 
may be more useful than others. 

C. STABILITY ON THE FREQUENCY PARAMETER PLANE 

Figure 3.8 shows the stability on the root parameter plane 
of the above system. The curve of ? = 0 is obtained from: 

1.248 - 0.4902o)^ + 5.541o3^ + 0.19968o)^ 

a = 5 

o)(25a)^ - 62.5) 

4 on 2 
o _ 0) - 30w 

P 2 

25o)^ - 62.5 

In the typical system, the region of stability on the 
frequency parameter plane can be determined by using the 
above relations. 

The stability boundary (?=0) curve has been added on 
Figure 3.1, through Figure 3.6 and the shaded areas are the 
stable region. Synthesis and design with the frequency para- 
meter plane are considered only for stable systems. 

Each Figure has a different interpretation. Figure 3.1 
plots a versus B with magnitude constant at -3 db , and frequency 
as the family parameter. By choosing an operating point in the 
stable region at the intersection of two frequency curves, one 
guarantees both stability and bandwidth. Figure 3.2 uses 
frequency as the constant parameter and magnitude as the 
family parameter. This permits optimum choice of a and B 
for desire magnitude and frequency. 
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The other curves can be considered in the same way. At 
any point on the Figures in the stable region, the parameters 
at that point can use to design the system and obtain values 
of magnitude and frequency as desired. 




Figure 3.1 Frequency Parameter Plane with Stability Limits. 

Plot of a vs 3 with M constant at -3 db and w 
as family parameter. 
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Figure 3.2a Frequency Parameter Plane with Stability Limits. 

Plot of a vs 3 with w constant at 10 rad/sec 
and M as family parameter. 
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Figure 3.2b Frequency Parameter Plane witJi Stability Limits. 

Plot of a vs S with co constant at 12 rad/sec 
and M as family parameter. 
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Figure 3.2c Frequency Parameter Plane with Stability Limits. 

Plot of a vs 3 with w constant at 14 rad/sec 
and M as family parameter. 
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Figure 3.3 Frequency Parameter Plane with Stability Limits. 

Plot of a vs M with w constant at 12 rad/sec 
and B as family parameter. 
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Figure 3.4a Frequency Parameter Plane with Stability Limits. 

Plot of a vs u with 3 constant at 3=4 and M as 
family parameter. 
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Figure 3.4b Frequency Parameter Plane with Stability Limits. 

Plot of a vs u with g constant at g=5 and M 
as family parameter 
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Figure 3.4c Frequency Parameter Plane with Stability Limits. 

Plot of a vs w with g constant at 3=6 and M 
as family parameter. 
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Figure 3.5 Frequency Parameter Plane with Stability Limits. 

Plot of a vs M with 3 constant at 3 = 5 at o) as 
family parameter. 
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Figure 3.6 Frequency Parameter Plane with Stability Limits. 

Plot of a vs w with M constant at -3 db and 3 
as the family parameter. 
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Figure 3.8 Root Parameter Plane Plot- 
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IV. 



nonlinear system techniques 



Generally, nonlinear systems consist of linear parts 
which are described by linear differential equations and 
nonlinear parts which are described by nonlinear functions. 

If the system Is nonlinear with an applied sinusoidal Input, 
the signal at various points In the system »ay or may not be 
sinusoidal. In general the signals will be periodic, but 
may contain sub-harmonics and/or harmonics. 

The most widely used method In the stability analysis 
of nonlinear control systems linearized by the describing 
function technique Is the Nyquist diagram. With the concept 
of the variable critical point, which depertfc upon the frequency 
and signal level at the Input of nonli nearf ty , the Nyquist 
diagram permits a simple stability analysis and Investigation 

of sustained oscillation. 

However, the Nyquist diagram Is Inconvenient to apply to 
multiloop control systems. Difficulties also arise with the 
adjustment of the system parameter, since t!he diagram should 
be plotted each time a parameter Is changed.. Thus, this 
method can hardly be applied to design problems In which the 
describing function does not appear as a variable gain factor. 

Other methods may be used, namely the IHurwItz criterion. 
Root-locus techniques. Bode diagrams, Nichoils charts, 

Mitrovic's method, approximation method. Popov method, and 
phase plane techniques, etc. which can be used to investigate 
the nonlinearity of the system. However the parameter plane 
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method is more convenient since the nonlinearity values can 
be represented by the parameters. 

Consider the sinusoidal system; 

D(jw,a,BJ 



Then 



M 



= |T(joj)| = 



N ( j (0 , g , g ) 
D ( j (jj , a , S ) 



Define one parameter as B to be constant and the other para- 
meter a represents the value of the nonlinear element. So; 

M 



N(.i>.g) 
D( jco,g) 



The equation has three variables: magnitude, frequency, and 
the nonlinear parameter a, which the parameter plane method 
is to plot. 

The block diagram of the nonlinear system may be the 
foil owi ng : 



N 




a is the nonlinear element. 

6(s) represented the linear system. 
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Nonlinear elements such as square law devices, coulomb 
friction, dead space, saturation, backlash, relays, etc. may 
be represented by their describing function. The parameter 
is then used to represent any adjustable parameter in the 
describing function, and the resulting plots on the frequency 
parameter plane indicate the effect of the nonlinearity on 
the system behavior. 
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V. DYNAMIC SYSTEM DESIGN 



A. INTRODUCTION 

In the design of compensators for linear or nonlinear 
systems, linear system design to make the system stable uses 
the root locus, Bede plot, Nyquist and Nichol chart. However 
the parameter plane method is more convenient for both linear 
and nonlinear systems. In classical methods of the root 
locus and frequency domain (Bode diagram), the compensators 
are introduced to reshape the corresponding graphs of the 
original system in order to achieve desired system response 
cha racteri sti cs . 

In the parameter plane approach the direct effects of 
the additional blocks and their connection Is to change the 
characteristic equation, so that the additional adjustable 
parameters are conveniently distributed amofig the coefficients 
of the characteristic equation. This enables us to have 
control over the characteristic roots by varying the introduced 
parameters . 

The parameter plane is used here to investigate and 
design the system for antiroll stabilization of a ship by 
means of passive tanks. The size and properties of the tank 
are designed. The problem is extended by Including the pump 
system thus converting the system to an active stabilization 
system . 
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B. THE DYNAMICS OF A SHIP IN ROLL MODE 

A mathematical description of the dynamic motion of a 
ship requires six degrees freedom, three equations of force 
and three moment equations, these six equations are too 
complex for analytical solution. If only the rolling of 
the ship is considered, it may be modeled as a pure rota- 
tional motion about the longitudinal axis of ship, with no 
moments about the Y and Z axes. 

The ship in rolling, is analogous to a pendulum. The 
dynamic equation of a pendulum is: 
m 4) + B 4) + 4) = F 

where m = mass of the pendulum. 

B = damping coefficient, 

g = gravity of the earth. 

1 = the length of the pendulum. 

<j) = the instantaneous angle of pendulum. 

F = the excite force. 

Including hydrodynamic damping, the dynamic rolling 
equation of a ship is the same as the pendulum. 

J^$+B (j) + K <})=Ky 

5 y S ^ S ^ Y 

where = total moment of inertia around axis of roll. 

Bg = hydrodynamic damping coefficient. 

Kg = static righting moment coefficient. 

Y = angle of drift, radians. 

= drift-angle coupling coefficient. 

<{) = instantaneous angle of rolling. 
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The equation can be written as: 



$ + 2?s ^ 



B. 



where t; = damping factor of the ship. = H — ^ 

V *^5 s 






Wg = natural frequency of the ship roll 
K 

When the ship is operating in a sea state, there is a > 
multitude of forcing frequencies present and the integrated 
response to all these frequencies must be considered. Further 
experiments are therefore carried out with the model at zero 
speed in irregular beam waves. The rolling of the ship is 
caused by the force due to the wave slope. Then the equation 
becomes : 

$ + 2 Cs “s ^ “s ” “s 
where the magnitude of wave slope. 



C. MODEL OF A SHIP 

The model of ship used for the example, to simulate and 
design the tank stabilizer is one given in Reference [6], the 
data and characteristics of the ship were the following; 
Displacement weight w = 936 tons. 

Metacentric height GM = 0.73 meter. 



Ship moment of inertia 



= 1290 tons-m-sec . 



Hydrodynamic damping coefficient = 60 tons-m-sec. 

Static righting moment coefficient = 684 tons-m. 
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The mass moment of inertia of the tanfe stabilizer (J^) 
increased the value of but it is very small and can be 
considered negligible when compared to 0^, for convenience 
use the value of as above. 

The beam of the ship can be found by. 



2it/GM _ 2ir/CT 




where C = a coefficient dependent upon the moment of roll 
and approximate is 0.797 sec/m^„ 

Therefore the beam of the ship can be cal cuiT ated and the 
result is: 

B = 9.25 meters. 

D. PASSIVE TANK STABILIZATION 
1 . Anal ysi s 

A passive tank can only be designie^d for optimum 
performance for one ship condition of di spTacement and meta- 
centric height. When the ship condition varies, the free 
surface correction y varies, together with ftihe quantity of 
water required. It is therefore necessary to check the 
performance of the stabilizer under differemt ship loading 
conditions and it could be that a compromise design has to be 
adopted which is acceptable for all conditions. 

The free surface loss is generally expressed as: 

i = 5 w B 

V 48 A h 
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where i = the second moment of inertia of the free surface 
7 = the volume of displacement of the ship, 
w = the weight of active fluid. 

A = displacement of the ship. 

B = beam of the ship, 
h = depth of fluid in tank. 

— = usual notation for free surface correction. 

St 

The tank size parameter is limited primarily by 
static stability requirements. However, the tank should be 
designed with as large a free surface as possible consistent 
with the requirements of static stability. Choosing values 
of y is equivalent to varying the tank size. 

Since the pendulum system is analogous to ship 
rolling, the effect on the ship of a hydrodynamic tank 
stabilizer is analogous to the addition of a second pendulum 
The combination of two pendulums in series, may be called a 
double-pendulum system as shown in Figure 5.1. 




Figure 5.1 Double-Pendulum System. 
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The first pendulum represents the ship rolling 
system, and the second pendulum represents the dynamic tank 
system. The length of the pendulums should be the same, so 
that the periods are nearly equal which means that both 
natural frequencies of the ship rolling and the tank are 
approximately the same. When the ship mass has been swung, 
the free oscillations are observed. When the tank mass has 
been attached, the masses of both ship and tank can be swung 
and the interplay and exchange of energy between the tank 
and the ship can be seen. 

The tank stabilization system uses water or some 
fluid mass device. There are many types of tanks. They can 
•be designed as Doughnut tanks, completely filled sea ducted 
tanks, partially filled sea ducted tanks, LE-tube tanks and 
free surface tanks. These tanks systems were suggested by 
Or. J. H. Chadwick Jr. Typical tanks such as free surface 
and U-tube tanks have been suggested by G. «J . Goodrich and 
Williams C. Webster respectively. Both types are described 
and designed in this chapter. 

A block diagram of ship and tank system is; 
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For a passive tank system the "Control block" 
does not exist, or may be considered a unit gain. Thus the 
block diagram for a passive system becomes: 



’I' 





SHIP 


<P 


J 

e 










TANK 









The relationship of the ship with a tank stabilizer 
can be represented as a double pendulum. The equations are: 



where 



Jgj + B^4> + Kg(|) + 

^ ®t® + = 0 

= the moment of inertia of the ship rolling. 

= hydrodynamic damping coefficient of ship. 

= the static righting moment coefficient of 
ship. 

= the moment of inertia of the tank. 

B^ = the damping coefficient of the tank. 

= the static equivalent coefficient of the 
tank . 

'^st ~ coupling moment of inertia of ship/tank. 

= the coupling static coefficient of ship/tank. 
F^ = the force insert in ship as the wave slope. 
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2 . Free Surface Tank Design 

The free surface tank as Figure 5.2 shown, is 
advantageous in decreasing the magnitude gain of rolling, 
and also decreases stability of the ship since the free 
surface loss is high. 



Figure 5.2 Free Surface Tank. 

The equations of the ship and tank can be rewritten 
as: for the ship 






for the tank 




where = natural frequency of the ship roll 



0 )^ = natural frequency of the tank. 




Cg = damping factor of the ship roll. 
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, = damping factor of the tank. 

2 

g = gravity of the earth = 9.8 m/sec. 

d = the distance between the water level in the 
tank and the roll center. 

ij; = wave slope function. 

p = the free surface correction factor. 

_ i _ 1 

■ V GM 2AGM 

Both equations have the negative sign in the coupling 
terms since the restoring forces are in opposite directions. 

If the distance d is positive which means the level of water 
in the tank is higher than the rolling center, the magnitude 
gain is less as shown in Figure 5.3, so for convenience set 
d equal zero. 

Generally, the period of rolling of the tank (or 
natural frequency) must be designed nearly equal to the 
period of the rolling of the ship. The simplified equations 
become : 

•• • 2 2 2 

<|) + a^(|) + 0)^ (j) - p e = 

0 + a^6 + = 0 

where = ^^s“s 

“t ' ^'t“t 

“n ' “s ' “t 

By manipulation and transform, the system trans- 
mission function of the rolling angle due to the wave slope 
force is: 
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. w^(S^+a.S+o)^) 

= D t n 

S^+(a^+ag)S^+(2w^+a^ag)S^+o)^(a^+ag )S+o)J( 1 -y) 



The solution of the equation can be expressed as 

an amplification factor A(u) and e the phase angle lag of 

the forcing function. 

2 



A (u) = 



0 ^^ 



/x^' 



+ Y‘ 



where X 



tan-1 Y 



( 2 2. 

= (o)^ - 0) } 



4/ 2 2s 
y w^(to^ - u ) 

~rz 2 \2 . 2 2 

(co^-to ) + a^w 



= o)[aL + 



yatcon 



-] 



s / 2 2s 2. 2 2 
(w -u ) +a.w 
11 u 

These solutions can be used to investigate the 
effect of variation of the tank parameters and y. Calcu- 
lations can be carried out using variation ©f each parameter 
in turn, keeping the other constant. Adjust both parameters 
to obtain the minimum magnitude gain. These values of para- 
meters are used to design the tank stabilizer. 

The computer simulation program used is due to 
Glavis [7]. The data of the ship in the previous section 
were substituted into the equation of ship and tank stabilizer, 
two unknown parameters in the equation are and y as variables 

Recal 1 

4(s) = T(s) = 






S^+(a^+a^)S^+( 2a)^+a^a^)S^+u^(a^+a^)S+a)^( 1 -y) 
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Numbering the data of the ship and parameters the 



equa ti on becomes : 



T(s) = 



0.53023S^+0.53023a^S+0.28115 



S^+(0.0465+a^)S^+(l . 06047+0 . 046 5a^) S ^+( 0 . 02466+ 



0.53023a^)S+0.28115(l-y) 
For the sinusoidal signal the equation becomes: 



2 . 



T(jo)) = -4 



0 . 281 1 5-0 . 53023oj +j0 . 53023a^co 



Then 

= 



(u’-(l .06047+0 . 04651aj)M+0.2Bl 15(1 -y)+j [(0.024662 
+0 .53023a.^)to-(0 .0465+a.j.)(o^] 

( 0 . 281 1 5-0. 53023w^)^+(0.53023a^(o)^ 

[o)^-(l . 06047+0.04651 a^) (0^+0. 281 15(l-y)]^+[( 0.024662 



+0.53023a^)w-(0.046 5+a^)(o'^] 



Since the size of the tank depends on the free 
surface coefficient (y) , if y is large, the tank will be big 
and the stability of the ship will be decreased because the 
metacentric hei ght decreases. So the magnitude gain of 
rolling in the region of the natural frequency decreases but 
at low frequency the magnitude gain increases. 

Consider the computer results. Figure 5.4 shows 
curves of the tank damping coefficient versus frequency, with 
magnitude as curve index and the free surface coefficient is 
constant. On Figure 5.4 the free surface coefficient is 0.05. 
If a (tank damping coefficient factor) is approximately 0.25 
the height of the resonance peak is 13 db , which is the 
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smallest value obtainable for these conditions. This occurs 
at a frequency of about 0.72 rad/sec. At a frequency of 
about 0.18 rad/sec, the magnitude is only 1 db, and is even 
smaller at low frequency. 

On Figure 5.5, the free surface coefficient is 0.1, 
the line of constant tank damping coefficient factor is 0.35 
to give a resonance peak of about 11 db. which is the minimum 
obtainable. Note, however, that at lower frequencies the 
magnitude is higher than for the condition of Figure 5.4. 

If the free surface is set equal to 0.2 as on 
Figure 5.6, the tank damping coefficient should be 0.45, 
which gives 9 db for the resonance peak and 2 db at low 
frequency (0.5 rad/sec). Thus, if the free surface coeffi- 
cient is increased, the magnitude in the region of low 
frequency also is increased, though the height of the reso- 
nance peak is decreased. 

The curves of Figure 5.7, Figure 5.8 and Figure 5.9 
give tank damping coefficient factor versus magnitude. 

Frequency is the curve index, and the free surface coefficients 
are 0.05, 0.1 and 0.2 respectively. These Figures correspond 
to Figure 5.4, Figure 5.5 and Figure 5.6. By inspection they 
predict the same optimum values for the damping coefficient 
factor. 

Figure 5.10, Figure 5.11 and Figure 5.12 give the 
Bode diagrams (plot of the magnitude versus frequency) with 
the tank damping coefficient factor as curve index and the 
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free surface coefficient factor is constant. In this case 
the Bode diagram cannot show the best value of the tank 
dampipg coefficient factor unless that value has been plotted. 
However the Bode diagram is good for consideration of the 
values of the free surface coefficient for which the curves 
show different magnitudes at all frequencies. 

By inspection of the curve families the best values 
of the tank damping coefficient factor and the free surface 
coefficient factor are: 

the tank damping coefficient factor = 0.25 

the free surface coefficient factor = 0.05 

and/or 

the tank damping coefficient factor = 0.35 

the free surface coefficient factor = 0.1 

Note that the magnitude gain value for low fre- 
quencies increases with increases of the free surface (y), 
which means increasing tank size. This results from the 
losses in roll stiffness due to the tank free surface effect. 
It is also obvious that the roll response at low frequencies 
increases as tank size increases. 

For the free surface tank (as reference [8]), the 
relation between the parameters and the tank size is; 

1 b^ 

^ 2A6M 

where 1 = the length of the tank. 

b = the width of the tank can be equal or less than 
the beam of the ship. 

A = the displacement of the ship in tons. 
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6M = metacentric height 

The other suggestion equation is: 



0 ) 



h = 



2 2 
t b 



2 „ 
ir g 



where h = the depth of the water in tank. 

If the width of the tank is assumed to be equal to 
the beam of ship, the depth of the water in the tank is equal 
to 0.5 meter. If the selected free surface coefficient factor 
is 0.1. 

1 = 0.82 meter. 

If the free surface coefficient factor is equal to 
0.05, the length of the tank is reduced to: 

1 = 0.41 meter 

which can save space in the ship. And the other parameters 
such as J^, B^, and can be calculated. 

The time response to an initial condition is shown 
in Figure 5.13, for y = 0.0, 0.05, 0.1, and 0.2, and = 0.0, 
0.25, 0.35 and 0.45 respectively. The curve for y = 0.0, 
without antiroll tank stabilizer is very oscillatory, and is 
damped only by the damping of the ship itself. The other 
curves include the tank stabilizer, so the damping is higher, 

■ 'and the curves reach steady state sooner. The time constant 
of the undamped system is about 42.5 seconds. With stabilizer 
added the setting time can be halved. 
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Figure 5.14 shows the time response to a sinusoidal 
input. The parameters have been varied as for Figure 5.13. 

The frequency of the input is set at about the natural fre- 
quency (0.728 rad/sec). For the ship without tank stabilizer 
the magnification is seen to be excessively large. When the 
tanks are added the system response is seen to be greatly 
reduced in magnitude. 

3. U-Tube Tank Design 

The U-Tude tank as designed in this chapter is 
assumed to be a rectangular tank. The mutual inertia term, 
Jg|.> which has been set to zero in preceding sections, is not 
necessarily zero. The value of depends on the accelera- 
tion of the water relative to the tank, and this depends on 
the location of the tank. If the cross-duct passes below 
the center of rotation, is a positive number which becomes 

smaller as the tank system is placed higher, as shown in 
Figure 5.15a and b. 

Consider that the two parameters and are 



given by: 



J 



t 



J 



St 



pAq 



ds 

0 



pA. S D. ds 

0 / 



where p = mass density of the fluid in tank. 

Aq = cross-sectional area of vertical leg of the tank. 
R = lever arm of the tank. 
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Dg = perpendicular distance from "virtual" center of 
rotation to the tangent at any point on the 
trajectary . 

Ag = the cross-sectional area at any point S. 

These two equations can be related as: 



Jst 

h 0 " ''o 



= 2h + 



2 

s s 



and h = the height of the water in the tank 



Note that can be negative if the tank is placed 
above the center of rotation. So can be zero if: 

-0 

s A. 



Recal 1 

■’st* + ''st* ^ ■’t® ^ = 0 



By using the definition as before the transformed 
equations become. 

(S^+a^S+w^)<J) + (- ^ S^- 3^)0 = 



(- S^-w^)(J) + (S^+a.S+to^)e = 0 

u ^ n L n 
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2 

Define = 3 (=yw, as in the free surface tank) 
'^s ^ 



St 



2 D. A. 



= X = 2h + 



Otj 



= a 



then 



St 



ex 

2 



0) 



n 



Using Cramer's rule, the ship roll transfer function 
due to the wave slope force is: 



|(s) 



w^(S^+atS+w^) 



( 1 -•^^)S^+(ag+a^)S'^ + ( 2w^+a^ag-2XB)S'^+w^(a^+ag)S+wJ^ 
0) 



2 . 2 , 



Inserting the values of the parameters of the ship, 
the equation becomes: 




0.530235^+0. 53023aS+0. 28115 

[(1 -1 .886X^3) S^+( 0.046 5+a)S^+( 1 .06047+0.0465-2X3)5^+ 

( 0. 02466+0. 53023a) S+0. 281 15-0. 530233 



then 

^2 ^ (0.281 1 5-0. 53023oj^)^+(0.53023aco)^ 

[(1-1 .886X^3)w^-(l . 06047+0. 0465 -2X3 )w^+0. 281 1 5-0. 530233]^ 

+ [(0. 024662+0. 53023a)w-(0.0465+a)w^]^ 



Figure 5.16a, b and c are plots for X=3 and 3=0.05, 
there are two resonance peaks in the magnitude curve, one 
for the ship itself and the other for the tank. The value of 
a at the minimum peak magnitude is about 0.45 and only one 
peak magnitude is obtained. The point P in Figure 5.16c 
corresponds to the line of -2 db in Figure 5.16a. 
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If the tank is placed as high as possible, the value 

of X may be zero or negative. If X=0, the result is the same 

2 

as for the free surface tank with 6=yw^. The plot of X=-1.0 
and 3=0.05 is shown as Figure 5.17a, b, and c. The points P 
and Q on Figure 5.17c are the intersection of the a lines, 
these points correspond to the lines of the frequencies 0.6 
and 0.9 rad/sec on Figure 5.17a. In this frequency region, 
if the value of a is decreased, the magnitude also decreases. 
Figure 5.17b shows that minimum peak value of magnitude is 
obtained with a=0.7. 

The U-Tube tank design equation is: 

Kj = 2CgAo 

where P = mass density of the water, 
g = gravity of the earth. 

Aq = cross-sectional area of each vertical leg of the 
tank . 

R = distance from ship center line to center of a tank 
vertical leg. 

The product of p and g is unity for the water, so both 

can be eliminated from the equation. 

And the cross-sectional area of the channel between 

tanks is given by: 

. 2 




Figure 5.3 Frequency Parameter Plane (Bode Diagram). Plot 
of the effect of the water level in Passive 
free-surface tank. 
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Figure 5.4 Frequency Parameter Plane for Free'^surface Tank. 

Plot of Tank Damping Factor vs Frequency with 
Free-surface , p=0.05 and M as curve index. 
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Figure 5.5 Frequency Parameter Plane for Free-^surf ace Tank. 

Plot of Tank Damping Factor vs Frequency with 
Free-surf ace , y=0.1 and M as curve index. 
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Figure 5.6 Frequency Parameter Plane for Free'^surface Tank. 

Plot of Tank Damping Factor vs Frequency with 
Free-surface , y=0.2 and M as curve index. 
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Figure 5.7 Frequency Parameter Plane for Free-surface Tank. 

Plot of Tank Damping Factor vs Magnitude with 
Free-surface, y=0.05 and w as curve index. 
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Figure 5.8 Frequency Parameter Plane for Free-surface Tank. 

Plot of Tank Damping Factor vs Magnitude with 
Free-surf ace , y=0.1 and u as curve index. 
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Figure 5.9 Frequency Parameter Plane for Free^surface Tank. 

Plot of Tank Damping Factor vs Magnitude with 
Free Surface, y = 0.2 and w as curve index. 
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Figure 5.10 Frequency Parameter Plane for Free-.surface Tank. 

Bode Diagram, Plot of Tank Damping Factor as 
curve Index with Free-Surf ace , y = 0.05. 
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Figure 5.11 Frequency Parameter Plane for Free-surface Tank. 

Bode Diagram, Plot of Tank Damping Factor as 
Curve Index and Free-surface , y = 0.1. 
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Figure 5.12 Frequency Parameter Plane for Free-surface Tank. 

Bode Diagram, Plot of Tank Damping Factor as 
Curve Index and Free-surface, p=0.2. 
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Figure 5.13 Initial Condition Response with Different Value 
of Tank Damping Factor and Free-surface. 
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Figure 5.14 Sinusoidal Response (co=.728 RAD/SEC) with Different V.alue of Tank 
Damping Factor and Free surface. 





(b) 

Figure 5.15 U-Tube Tank 
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Figure 5.16a Frequency Parameter Plane for U-Tube Tank. 

Plot of Tank Damping Factor vs Frequency 
with X = 3, 3=0.05 and M as Curve Index. 
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Figure 5.16b Frequency Parameter Plane for U-Tube Tank. 

Plot of Tank Damping Factor vs Magnitude 
with X = 3, 3=0.05 and w as Curve Index. 
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Figure 5.16c Frequency Parameter Plane for U-Tube Tank. 

Bode Diagram, Plot of Tank Damping Factor 
as Curve Index with X = 3 and 3=0.05. 
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Figure 5.17a Frequency Parameter Plane for U-Tube Tank. 

Plot of Tank Damping Factor vs Frequency 
with X = -1 , 3=0.05 and M as curve index. 
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Figure 5.17b Frequency Parameter Plane for U-Tube Tank. 

Plot of Tank Damping Factor vs Magnitude 
with X = -1 , 3=0.05 and u as Curve Index. 
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Figure 5.17c Frequency Parameter Plane for U-Tube Tank. 

Bode Diagram, Plot of Tank Damping Factor 
as Curve Index with X = -1 and B=0.05. 
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This is for the rectangular U-Tube tanks. With the 
equations as above, the values of cross-sectional area of 
the vertical leg of the tank, Aq, the cross-sectional area 
of the channel, , can be calculated and tank designed. 

The depth of the water in tank, h, and the distance 
between the rolling center and the channel of the tank, , 
can be calculated by: 



In the case of X = 3, if the tank is located below 
rolling center, h should be less than 1.5. 

E.g., if X = 3, g=0.05, the tank damping factor is 
0.45. Then the tank size is: 



if R 



Aq R- 



= 4, 



•^t 4 

^ 2Pg ^ 2Pg ' 3^.25 m 

the width of the tank is about 1 meter. 



The ref ore 



Aq = 2.013 m^ 

2 

and if h = 1.3 meters, the result of A„ is 0.4686 meter 

s 

and = 0.859 meter, 
s 



E. ACTIVE TANK STABILIZATION 
1 . Analysis 

Acti vated-anti rol 1 tanks are considered in detail. 
Particular attention is paid to the design of the tank control 
system. The equations of motion are the same as those of the 
ship with passive tank stabilizer, but must be modified to 
include the pump used in the activated tank. Active tank 
control can best be selected on the basis of minimizing the 
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response of the ship to an impulsive roll velocity, and it 
will be shown that the active tank is significantly better 
in reducing the roll of ship than the ordinary passive antiroll 
tank. The active tank stabilizer studied uses the U-Tube tank 

f 

design. 

Active tank analysis parallels that of passive tanks 
and/or the double pendulum. The tank system is activated by 
inserting a force (such as a pump system) into the second 
equation of passive tank system. The equation can then be 
written as: 

ship: Jg(}) + Bg({) + Kg({) + J^^e’ + 

tank: + J^e + B^e + K^e = Fp 

where Fp = the force due to pump system. 

The effect of the roll sensing control on the force 
produced by the pump system can be expressed as: 

K„(6n+6,S+6.jS^)({) 

F (s) = P ^ ^ — 

1+ tS 

where Kp = the gain of the control pump. 

6-| = roll angle sensitivity. 

62 = roll rate sensitivity. 

6^ = roll acceleration sensitivity. 

T = time delay constant. 

For convenience, we set x =0 in the computation, since 
the time delay is small. 
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The active tank system can be shown in block diagram 



form as : 




where F^(s) = transfer function of the ship rolling due to 

wave slope alone which is the same as for the 
passive tank system. 

P^(s) = transfer function of the ship rolling due to 
pump alone, which can be expressed as; 



(l-^)S^+(a^+a^)S^+(2a3^+a^a2-2XB)s2+(o^(a^+a^)S+(oJ-Bw^ 

F„ 



2 . System Design 

From the analysis in this section, the system per- 
formance can be studied using computer simulation. The pump 

• ' 

coefficient factor (K^ = is one variable parameter 

p p VJ^ 

which may be called a. It is varied as the pitch ratio of 
the propeller of the pump, and can be assumed to be varied 
from zero to one and zero to minus one. For convenience the 
actual system may be designed with only positive variation and 
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the direction of the water motion between tanks is controlled 
by valves. 



From the block diagram the system transmission 
function is: 



(l-^^^)S^+{a^+a^)S^+(2o)^+a^a^-2X6-K GjB)S^+(a)^aj.+ti)^a, 



f(^) = 






' 4^2 



-KpG2)S+a3n-3o3n-KpGiB 



where ^ Kp( G^ +G2S+G2S^ ) <)) 

Define X. =3 

3 = 0.05 

a. =0.4 



and 



K. 



= a 



These values for a passive tank gave a peak resonance 
of minimum magnitude at about 10 db at 0.6 rad/sec and is 1 db 
at the lower frequencies. Using these values the equation 
becomes : 



(s) = 



0. 530235^+0. 21 2092S+0. 281 15 



0. 1 51 3S^+0. 44655^+ (0. 7791 -0.05aG3)S^+(0. 237-0. 05aG2)S 



+0.25464-0 .05aG 



1 



Figure 5.18 shows the Bode diagram for G^=l, 
82=63=0. If a is more positive, the frequency and the magni- 
tude of the resonance peak are slightly decreased. At the 
low frequencies, the magnitude is increased. If a is more 
negative, the re suit is reversed. To obtain the mini mum 
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Figure 5.18 Frequency Parameter Plane for Activated Tank. 

Bode Diagram, Plot of Pump Coefficient Factor 
as Curve Index (G^=l, 62=62=0) and X=3. 
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magnitude gain, at the lower frequencies from 0 to 0.6 rad/sec 
a is -1, for the higher frequencies a must be +1. So, the 
pitch ratio of the propeller of the pump may be kep constant 
and the direction of the water motion reversed by switching 
the valves. Then the maximum magnitude is 7 db at 0.6 rad/sec 
and -1 db at the lower frequencies. 

The characteristics of the active device (pump, 
valves, sensor) can be treated as the variation of the value 
of the parameter a versus frequency. This value is required 
to switch at the resonant frequency as shown in Figure 5.19. 



1 

a 



0 



0.6 



0 ) rad/sec 



-1 



Figure 5.19 The variation of control device. 

(X = 3 , G^ = l , 02 = 63 = 0) 

For G-|=0.5, 03 = 1.25 and 03 = 2.4 the results are shown 
as Figure 5.20a, b and c. Figure 5.20a is a plot of a versus 
frequency, and magnitude is the curve index. When a is more 
negative the magnitude is seen to be smaller. If a=l the 
maximum peak magnitude is obtained at the frequency of about 
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Figure 5.20a Frequency Parameter Plane for Activated Tank. 

Plot of Pump Coefficient Factor vs Frequency 
with M as Curve Index (G-.=0.5, G^ = 1.25, 
63 = 2 . 4 ) and X=3. ' 
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Figure 5.20b Frequency Parameter Plane for Activated Tank. 

Plot of Pump Coefficient Factor vs Magnitude 
with 0 ) as Curve Index (G,=0.5, G., = 1.25, 
63 = 2 . 4 ) and X=3. ' ^ 
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Figure 5.20c Frequency Parameter Plane for Activated Tank, 
Bode Diagram, Plot of Pump Coefficient Factor 
as Curve Index (G,=0.5, G^=1.25, G^=2,4) and 
X = 3. ' 
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0.6 rad/sec. These results correspond to Figure 5.20b which 
shows the frequency curve. At the same frequency of 0.6 
rad/sec, both figures show maximum magnitude for all values 
of a from -1 to +1. Figure 5.20c is the Bode diagram and 
shows that there is no phase lag between ship and active tank, 
and if a is more negative the magnitude is smaller. So, at 
a=-l , it produces the smallest magnitude at all frequencies, 
with the peak equal 5 db at 0.6 rad/sec. 

This result seems to be the best since all the 
parameters of the pump system (pitch, sensing gain control, 
and valves) are fixed. 

From both results as mentioned above, the system 
can be designed as a linear device. 

If there is no limit on the height at which the 
tank can be located, it is possible to choose the height, so 
that X= 0 . If this is done, by inspection the value of 3=0.05 
and a^=0.3 are obtained. Substitute these values into the 
equation of the system which becomes; 

/..X _ 0.53023S^+0.1 59069S + 0. 28115 

(s) = -j = P 

S^+0.3465S'^+(1 .07436-0. 05K G, ( S^+( 0 . 18373-0 . 05K G,)S 

p 3 P 4 

+ 0. 25464-0. 05KpG.| 

Figure 5.21a and b show the results with G-|=l and 
62 = 62 = 0 . Figure 5.21a is a plot with magnitude as curve 
index; point Q is the maximum magnitude with a=-r0.65 at a 
frequency of 0.728 rad/sec. If a=0.15, the maximum magnitude 
of 9 db is obtained which is the minimum peak. Figure 5.21b 
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Figure 5 . 21 a Frequency Parameter Plane for Activated Tank. 

Plot of Pump Coefficient Factor vs Frequency 
with M as Curve Index (G-|=l, 62 = 62=0) and X= 0 . 
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Figure 5 . 21 b Frequency Parameter Plane for Activated Tank. 

Bode Diagram, Plot of Pump Coefficient Factor 
as Curve Index (G.|=l, 62 = 62=0) and X= 0 . 
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shows the Bode diagram. To obtain the minimum magnitude, in 
the frequency region between 0.6 and 0.85 rad/sec the value 
of a must be maintained at +1 and for the other frequencies, 
at -1. There are two magnitude peaks of approximately 7 db 
at the intersection of a=+l and a=-l. For the lower frequen- 
cies the magnitude is -1 db. 

From the above result the active device parameter 
value, a, can be plotted versus frequency as in Figure 5.22. 
This figure shows two values of a at the frequencies between 
0.6 and 0.85, the straight line corresponds to the minimum 
magnitude, and the curved line corresponds to a constant 
magnitude in that region. 




Figure 5.22 The Variation of Active Device 
(X=0, G^=l, 62 = 63 = 0 ) 

Figure 5.23 shows the Bode diagram of the system 
for 6-|=0.5, 63 = 1 . 25 , and 63 = 2,4. For minimum magnitude with 
the control device, the required variation in a is shown in 
Figure 5.24. 
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Figure 5.23 Frequency Parameter Plane for Activated Tank. 

Bode Diagram, Plot of Pump Coefficient Factor 
as Curve Index (G.,=0.5, G« = 1.25, G^=2.4) and 
X=0. ' 
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Figure 5.24 The Variation of Control Device 
(X=0, G^=0.5, 02=1.25 23 = 2 . 4 ) 

Figure 5.25, which is the Bode diagram, compares 
the results for an undamped ship and for a ship damped by a 
tank. All have only one resonance peak. The magnitude of 
roll in the undamped ship at resonance is the highest, damping 
clearly reduces the amplitude. However, in the passive tank 
system the magnitude at lower frequencies is slightly higher, 
and the curves show that the U-Tube tank is more advantageous 
than the free-surface tank. The active tank system gave the 
most advantageous behavior as shown. 
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Figure 5.25 Compare the Results in Magnitude Gain. 
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VI. CONCLUSION 



Analysis techniques can be employed effectively in system 
design particularly when an exact analytical solution is not 
available and where a trial and error method is the only 
solution possible. 

The frequency parameter plane techniques have been used 
in this thesis to analyze and design a system. The system 
elements which may be linear or nonlinear were represented 
by adjustable parameters. To select a stable system, the root 
parameter plane was used and further mapped onto the frequency 
parameter plane. The goal was to obtain the best set of para- 
meters for a stable system. 

Rather than adjust parameters for every trial, one graph 
was made and the parameter value corresponding to the best 
curve was chosen. This appreciably reduced the number of 
trials and the computation time in the overall design. 

In the system which is presented the tank stabilizer 
design might not be optimal, however further effort on para- 
meter adjustments by iteration of the described technique 
would improve the design. 

The tank stabilizers of a ship were designed in this 
thesis by using the frequency parameter plane technique, they 
were very successful and effective. If the relationship 
between the designed parameters and tank size are true, the 
result of designing by using the method would be feasible and 
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believable. In passive tank design, the tank size depends 
on the free surface factor. The tank damping factor depends 
on the skin of the material of the tank, which would be tested 
to obtain its value from the analysis. The value of X in 
U-Tube tank would have to be chosen within the range which 
maintains ship stability. 

In the activated tank design, the gain control constants 
were varied by trial and error to adjust the damping. Phys- 
ically only one parameter can be varied, the pitch radio of 
the propeller of the pump. 
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